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Abstract 

We consider a coarse-grained model in which polymers under good-solvent conditions are repre- 
sented by soft spheres whose radii, which should be identified with the polymer radii of gyrations, 
are allowed to fluctuate. The corresponding pair potential depends on the sphere radii. This model 
is a single-sphere version of the one proposed in Vettorel et al, Soft Matter 6, 2282 (2010), and it is 
sufficiently simple to allow us to determine all potentials accurately from full-monomer simulations 
of two isolated polymers (zero-density potentials). We find that in the dilute regime (which is the 
expected validity range of single-sphere coarse-grained models based on zero-density potentials) 
this model correctly reproduces the density dependence of the radius of gyration. However, for the 
thermodynamics and the intermolecular structure, the model is largely equivalent to the simpler 
one in which the sphere radii are fixed to the average value of the radius of gyration and radii- 
independent potentials are used: for the thermodynamics there is no advantage in considering a 
fluctuating sphere size. 

PACS numbers: 61.25.he, 82.35.Lr 
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I. INTRODUCTION 



Polymer solutions are very interesting soft-matter systems, showing a wide variety of 
behaviors, depending on density, temperature, architecture, etc." Due to the large number 
of atoms belonging to a single macromolecule, simulations of polymer systems are quite 
challenging. For this reason, during the years many attempts have been made to develop 
coarse-grained (CG) models in which only some relevant degrees of freedom are retained, 
in such a way to reproduce some large-scale structural properties and the thermodynamic 
behavior.- In the simplest approach one maps polymer chains onto point particles interacting 
by means of the pairwise potential of mean force between the centers of mass of two isolated 
polymers.-"- Since the potentials are computed in the limit of zero polymer density, this ap- 
proach is limited to the dilute regime, in which many-body interactions^ii^ can be neglected. 
This limitation was overcome ten years ago,— by introducing pair potentials depending on 
the polymer density, thus allowing the model to reproduce exactly the thermodynamics at 
any given density. This work has paved the way to the use of soft effective particles to rep- 
resent polymer coils in complex situations such as in modelling colloid-polymer mixtures.— 
However, deriving density-dependent potentials requires full-monomer simulations at finite 
polymer density, which is what one would like to avoid by using CG models. Moreover, care 
is needed to derive the correct thermodynamics^"— and to compute free energies and phase 
diagrams.- 

These limitations can be overcome by switching to a model at a lower level of coarse 
graining, i.e. by mapping a long linear polymer to a short linear chain of soft effective 
blobs.—"— Quite recently, extending an older phenomenological appro ach,—'^ Ref. 2l| pro- 
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3, 



the potentials do 



posed a multiblob model. As in the models proposed in Refs. 
not depend on polymer density and are fixed by using structural data obtained in the limit 
of zero polymer density. However, at variance with the other approaches, the radius of each 
blob, which should be identified with its radius of gyration, is not fixed but is allowed to 
fluctuate. The idea is very appealing and the model is, in principle, more realistic, since 
blobs are compressible: when the density increases, the average blob radius of gyration de- 
creases as it does in the original polymer model. Therefore, this model should provide a more 
accurate description of the structural properties of the polymers. However, its practical im- 
plementation is by no means straightforward. A phenomenological approach was proposed 
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in Ref. |2l|, in which intermolecular and intramolecular potentials appropriate for polymers 
under good-solvent conditions were determined by combining exact results for ideal chains 
and some approximate predictions for good-solvent polymers that allowed the inclusion of 
the local self-repulsion. 



In this paper, we wish to test the approach of Ref. |21| for polymers under good-solvent 
conditions. We consider a single-blob system and carefully compare the results obtained 
for the model with radii-dependent potentials with those obtained for the simpler model 
in which the blob size is fixed. The simplicity of the CG single-blob model allows us to 
compute the pair potentials from full-monomer simulations of two isolated chains, avoid- 
ing any approximation, thus allowing us to distinguish between the merits/demerits of the 
method from those of the approximations which are needed to implement it. Moreover, we 
can discuss two different approaches to the coarse graining: in the first one each polymer is 
represented by a soft "compressible" sphere located in the polymer center of mass (this is the 
most common approach when dealing with linear polymers); in the second one the position 
of the sphere coincides with that of the central monomer (star-polymer studies favor this 
second option). 

The paper is organized as follows. In Sec. [Ill we define the models we investigate. In 
Sec. mil we present our results: in Sec. IIII Al we discuss the third virial coefficient and the 
zero-density three-body forces, while in Sec. IIII Bl we compare the results for the compress- 
ibility factor, the center-of-mass (or polymer midpoint) pair distribution function, and the 
distribution of the radius of gyration in the semidilute regime. Finally, in Sec. [IV]we present 
our conclusions. 



II. THE MODELS 

In this paper polymers are represented as "compressible" soft spheres of radius a: two 
spheres of radii o"i and (72, respectively, interact by means of the pair potential V{ai, 0-2; 6), 
where h is the relative distance. The radius of each sphere is allowed to fluctuate. We assume 
that the normalized radius distribution for an isolated sphere is given by the function -P(cr). 
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For a system of L spheres in a volume V, we consider the partition function 



Z= dcriP{ai) . . .d(rLP{crL)QL{cri, . . . ,crL) 



Qiiai, . . . ^ai) = / d^bi . . .d^bi exp 
Jv 



-l3^V{ai,aj;bi 



i>j 



where 6j is the position of the z-th sphere and bij = \bi — bj\. Equivalently, as in Ref. 
can define a potential pVi{a) = — In P(o") and write the partition function as 
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daid^bi . . . daid^bi exp 



i>j 



(1) 



we 



(2) 



The distribution P(cr) as well as the potentials V{ai, o"2; b) are determined in such a way to 
reproduce the thermodynamics and the intermolecular spatial distribution of the polymers 
in the limit of zero polymer density. More specifically, we consider a polymer model, in 
which each polymer consists of monomers located in ri, . . .r^. As usual we define the 
radius of gyration of the chain as 



-Ti 



2N^ 



(3) 



and its average over all polymer configurations as Rg{N) = (rp. We indicate zero-density 
averages with a hat, so that Rg{N) is the average radius of gyration of a single isolated 
polymer. It is a function of and for A^ — )■ oo (scaling limit) it scales as 



Rg{N) =aN''{l + bN-^ + 



(4) 



where a and b are model-dependent constants and u and A are universal exponents. For 
good-solvent polymers they are known quite precisely:— u = 0.587597(7), A = 0.528(12). 
To define -P(ct) we consider the distribution of the radius of gyration of an isolated polymer. 



Q{s-N) = {S{rg-s)), 



(5) 



where is the statistical average over all the conformations of an isolated polymer made 
of A^ monomers. In the scaling limit A^ — > oo, the adimensional quantity Rg{N)Q{s] N) is a 
universal function of a = s/Rg, which we identify with -P(o"): 



Rg{N)Q{s;N) = P{a) + 0{N-^) 



(6) 
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Scaling corrections decay as ^, where A is the same universal correction-to-scaling ex- 
ponent appearing in Eq. (jl]). The function P(cr) satisfies 

daP{a) = l, [daa^P{a) = l, (7) 



the second equation being a direct consequence of the definition of Rg. To define the poten- 
tial, we first consider the model in which the CG spheres are located in the centers of mass 
of the polymers. In this case we first consider 

/3v{si,S2;r;N) = - In ^ '-, 8 

Q{si;N)Q{s2;N) 

where the average (■)o,r is over all isolated polymer pairs made of monomers, such that 
their centers of mass are in the origin and in r, respectively, i and 2 are the gyration 
radii of the two polymers, and f/intcrm is the intermolecular interaction energy. In the scaling 
limit — 7- 00, v{si, S2',r; N) converges to a universal function, which we identify with the 
pair potential of the CG model: 

v{su S2; r; N) = V{ai, ^2; b) + ©(AT"^), (9) 

where cxi = Si/Rg, 02 = S2/-R3, and h = r/Rg. Because of these definitions, the polymer 
center-of-mass pair distribution function in the limit of zero polymer density is the same as 
the zero-density pair distribution function between the centers of the spheres. Hence, the 
polymer thermodynamics is correctly reproduced in the low-density limit by the effective 
model. 

It is important to stress that the quantities -P(cr) and V{ai,a2]b) are universal, i.e. the 
same result is obtained by using any model that is appropriate to describe polymers under 
good-solvent conditions. One could use the well-known lattice self-avoiding walk model 



or any off-lattice mode 
bead-rod model of Ref. 
(see Ref. 27 
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appropriate to describe good-solvent polymers, for instance the 
We will use the lattice Domb- Joyce model with w = 0.505838 



for details on the model). In the scaling limit A^ — ?■ 00 such a model describes 
polymers under good-solvent conditions. Morever, for our particular choice of the parameter 
w, the scaling corrections proportional to A^~^ that appear in Eqs. ([6]) and are very 
small, so that we can obtain the asymptotic (scaling) functions ^(0") and V{ai,cr2;b) from 
simulations of chains of moderate length. The results we shall present are obtained by 
using chains with A^ = 600. We have also performed simulations with A^ = 2400: the 



b 

FIG. 1: On top we report the center-of-mass potentials pV{ai, a2', b): on the left we show a three- 
dimensional plot in terms of a = ai = (T2 and b = r/Rg, on the right we show the potentials for 
cji = 1.024 and several values of (T2, as a function of b. On bottom we show the same plots for the 
midpoint potentials Va./p(cji, C72; 6). 

corresponding results are fully compatible with those obtained using = 600, indicating 
the absence of relevant finite-length effects. In practice, we simulate two independent Domb- 
Joyce chains using the highly-efficient pivot algorithni^I (y^Q perform 1.25 x 10^ pivot trial 
moves on each of them) and determine numerically Q{s, N = 600) and v{si, S2] r; N = 600). 
Then, we define P((t) = Rg{N = 600)g(s, = 600) and V{ai, as; b) = v{si, S2; r; N = 600), 
with s = aRg{N = 600) (analogous relations hold for si and S2) and r = bRg{N = 600). 
Plots of the potential for several values of cxi and a2 are reported in Fig. [1] They have 
an approximately Gaussian behavior with V{o'i,a2;b = 0) increasing as the radii decrease. 
This is of course expected, since the smaller o"i and o"2 are, the more compact the two walks 
become. Hence repulsion should increase. For the same reasons, the range of the potentials 
decreases as the radii decrease. 
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In Ref. |2l|, on the basis of a heuristic argument, it was suggested that the pair potentials 



could be approximately written as 

VvBKicTi, aa; b) = e{af + ajy^^^ exp (-a ^ \ , (10) 

where a = 3/2 and e is a constant independent of oi and Oi- We find that this expres- 
sion works reasonably well when taking e ~ 4-5. To obtain an optimal approximation we 
determine a and e in such a way to minimize the functional 

/>00 /-OO /-fcrnax „ 

^{a,e)= da^ da2 d6 [P((Ti)P(a2)62 (e-'^^^'^^-'^^''') - e-'^^^^^^'^^-'^^''))] , (11) 

Jo Jo Jo 

where V{cri,a2',b) and -P(cr) are the quantities computed from full-monomer simulations. 
Since the Monte Carlo estimates of V for b > 3 are noisy, we excluded these values from 
the b integration, taking 6max = 3. The functional ffTTl) has been chosen on the basis of the 
expression of the second-virial universal combination^^ 

A2 = 2tt j daida2dbP{ai)P{a2)b^ (l - e^'^^^"^'"^'^)) . (12) 

We obtain 

a = 1.42, e = 4.42. (13) 

Note that the optimal value for a is quite close to the value a = 3/2 proposed in Ref. 
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In 



Fig. [2]we plot {af + (7|)^/^V^(ai, 0^2; b) vs R = b{af + o\)~^l'^ for several values of Ox and 02- 
The results are then compared with the phenomenological expression e e~"^^ , obtained using 
potential (fTOj) . We observe reasonable agreement for i? > 1, while significant discrepancies 
are observed for i? — )■ 0. However, this is exactly the region which does not contribute 
significantly to A2, hence to the thermodynamics. 

Up to now we have characterized the polymer position by using its center of mass: in 
definition ([9]) b represents the distance between the two centers of mass, expressed in units 
oi Rg-. b = r/Rg. This is the usual choice for CG models of linear polymers. On the other 
hand, when considering star polymers, it is much more common to consider CG models in 
which the star position is identified with the position of the ramification point, the center of 
the star.-i^"— If we view linear polymers as two-arm star polymers, it is natural to define the 
CG model using the central monomer as polymer position. In an exact mapping this choice 
would be ininfluential for the thermodynamics.- However, since all n-polymer interactions 
with n > 3 are neglected in the CG model, different choices of the point associated with the 




FIG. 2: Rescaled potentials {af + a2)^^'^V{ai,a2; b) as a function of R = b{a\ + a"^) ^^"^ for several 
values of cJi and a2- We also plot the function ee~"^^, with e = 4.42, q = 1.42 (VBK). 

polymer position give rise to CG models with different thermodynamic behavior. It makes 
therefore sense to study the different possibilities, with the purpose of understanding which 
is the optimal one. 

We thus also study two different CG models in which the polymer position is given by the 
central monomer rjv/2- We consider a model in which polymers are represented by identical 
spheres interacting by means of the potential of mean force 

Vmp(&) = -ln(e-^^—)o,r, (14) 

where now the average is over all isolated polymer pairs such that the central monomers 
(midpoints, MP) are in the origin and in r, respectively, and h = r/Rg. We also consider the 
model with compressible soft polymers interacting by means of the potential Vmp{ci, cr2] b) 
defined as in Eq. ([9]), where now b is the distance between the central monomers expressed 
in units of Rg. 

The potential Va/p(6) has been discussed at length in the context of star polymers. For 
6 — i- it diverges logarithmically as^ ln(l/6). An explicit parametrization has been given 



in Ref . 



35 



(see their results for a two-arm star polymer): 



VMp{b) = - In 
r 



(^)''e-^''^+exp(r7e-^^^)' 



(15) 



where 

a = 1.869, /3 = 0.815, 7 = 0.372, S = 0.405, r = 4.5. (16) 
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This parametrization is quite precise. For instance, we obtain 5.51 for the second-virial 
combination A2, which is very close to the polymer result^ A2 = 5.500(3). 

In the bottom panels of Fig. [T] we report the potential Vmp(o"i, o"2; &). It is interesting 
to observe that they all diverge logarithmically as 6 — ?■ 0, apparently with the same type of 
logarithmic behavior, VMp{cri,o'2;b) ~ 0.821n(l/6), for all values of di and (T2. 



III. COMPARISON OF THE MODELS 



In this section we compare the thermodynamic behavior of the models we have introduced 
in the previous section: 

a) the standard CG model in which polymers are identical soft spheres interacting with a 
potential which only depends on distance. We consider the case in which the polymer 
position is given by the position of the center of mass (model Mia) or of the central 
monomer (model Mlb). In the first case we use the accurate expression of the pair 



potential^^ given in Ref. |37|, in the second one we use Eq. (fTSjl with parameters (fT6|) . 



b) We consider the "compressible" CG model in which we use either the center of mass 
(model M2a) or the central monomer (model M2b) as polymer position. We also 
consider the model with pair potential ffTOj) with parameters given by (1131) (model 
M2c). 

The results will be compared with full-monomer results and, for comparison, with those 
obtained in the tetramer model (results will be labelled with "t") introduced recently in 



Ref. 
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A. Three-body interactions at zero density 

If we expand the compressibility factor in powers of the concentration c = L/V as 

Z = = 1 + B2C + B,c' + 0(c3), (17) 

kbJ- c 

the quantity A2 = B2/Rg is universal. An accurate estimate is^i A2 = 5.500(3). Since we 
have matched the center-of-mass or polymer-midpoint distribution function to determine 
the pair potential, all models should give the correct estimate of the combination A2. In 
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TABLE I: Virial-coefficient universal combinations for the models introduced in Sec. IIIII and for 
the tetramer model (t) of Ref. 



22l . We also report the universal asymptotic values for polymers 



(P)- 
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p Mia 


Mlb 


M2a M2b 


M2c 


t 




5.500(3) 5.4926(1) 


5.5109(1) 


5.5102(3) 5.5085(2) 


5.5738(3) 


5.597(1) 


A3 


9.80(2) 7.844(6) 


4.925(4) 


7.42(2) 4.43(2) 


7.22(2) 


9.99(2) 


A' 


10.64 7.844(6) 


4.925(4) 


8.015(5) 5.012(2) 


7.809(5) 


10.57(2) 




-0.84 





-0.59(2) -0.58(2) 


-0.58(2) 


-0.581(5) 



Table U we report the results for the models we consider. Differences are small, and are 
representative of the level of precision with which the models reproduce the polymer center- 
of-mass or midpoint distribution function. Note that the estimate corresponding to model 
M2c is very close to the correct one, indicating that expression ( ITOl) parametrizes quite 
accurately the Rg dependence of the potentials. 

Much more interesting is the comparison of the third virial coefficient, since it provides 
an indication of the accuracy with which the CG models reproduce the polymer thermody- 
namics in the dilute regime and also of the importance of the three-body forces which have 



been neglected. The universal combination = B^/R^ was computed in Ref. 



271 finding 



A3 = 9.80(2). (18) 

In order to determine A3, two contributions had to be computed. One contribution is the 
standard one, which is the only one present in monoatomic fluids and in fluids of rigid 
molecules, A'^ ~ 10.64, while the second one is a flexibility contribution ^3,/^ ~ —0.84 (it 



corresponds to —TiR~^ in the notations of Ref. 1271 ). The combination ^3 as well as the two 
contributions A'^ and A3J1 are universal, hence it makes sense to compare them with the 
corresponding results in the CG models. 

We have estimated ^3 for all models. The results are reported in Table [H Note that, 
while A3J1 vanishes in models Mia and Mlb, a nonvanishing contribution with the correct 
sign is obtained for the models with a fluctuating radius. On the other hand, the estimates 
of A3 for models with radii-dependent potentials are essentially equivalent to those in which 
the sphere radii are fixed: the estimates corresponding to models M2a and M2b are close 
to those of models Mia and Mlb, respectively. To be precise, discrepancies increase by 
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FIG. 3: Three-body potential of mean force /3V3(ri2, ri3, r23) for ri2 = ri3 = r23 = r, as a function 
of b = r/Rg. On the left we report results for models M2a, M2c, for the tetramer model (t) of 



Ref. 122, and the predictions of full- monomer simulations for the quantity associated with the center 
of mass (FMa); on the right we report the results for model M2b and the prediction of full- monomer 
simulations for the analogous quantity associated with the polymer midpoint (FMb). 

allowing the radii to fluctuate: the difference between the M2a (or M2b) estimate of ^3 and 
the asymptotic polymer result is larger than the discrepancy observed for model Mia (or 
Mlb, respectively). We can also compare the results for models M2a and M2c: the difference 
between the two estimates of is small, confirming that expression (fTOj) is reasonably 
accurate. Finally, the results show that the CG models in which the center of mass is taken 
as reference point are more accurate than those in which the central monomer is considered. 
In the latter case, ^3 is underestimated by approximately a factor of two. 

As a further check we compute the effective three-body potential of mean force defined 
byiS^ 

here f/^ is the intermolecular potential energy between molecules i and j and the av- 
erage (■)bi2,bi3,b23 is over the radii distributions of triplets of isolated spheres such that 
bjj = bj — hj, where hi is the position of sphere i. In models Mia and Mlb, the sphere 
radius does not fluctuate, so that (e-^^"-')b,^ = e'^^^^'^) and (e~^^i2-/3C/i3-/3C/23^j^^^ = 
e-/3i/(f.i2)-/3y(fei3)-/3V(fei3)_ It follows that 1^3(612, &13, ^23) = 0. For models M2a, M2b, and M2c 
the average is over the distribution P{a) and Uij should be identified with V{ai, af, bij). Po- 
tential ( fT9l) should be compared with the analogous polymer quantity in which (■)bi2,bi3,b23 is 
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the average over all conformations of triplets of isolated polymers such that 6j is the position 
(in units of Rg) of the center of mass of polymer i (this is the case relevant for models Mia, 
M2a, and M2c) or the position of the polymer midpoint (the relevant potential for models 
Mlb and M2b). The polymer three-body potential of mean force is universal in the scaling 
limit, i.e. it is model independent. 

We computed /3V3(bi2, bis, b23) for equilateral triangular configurations such that 612 = 
^13 = = b for models M2a and M2b. The results are reported in Fig. [3] and compared 
with the corresponding quantities (FMa and FMb) computed in full-monomer simulations. 
They were obtained by considering triplets of Domb- Joyce walks made of = 600 beads. 
We used the pivot algorithm and performed 2.5 x 10^ pivot trial moves on each of them. We 
also performed simulations with = 2400, verifying the absence of finite- length effects. 

For model M2a, we find PV^lb, b,b) = for 6 > 1.2: for these values of b, model M2a is 
not different from the simpler model Mia. In particular, it does not reproduce the repulsive 
maximum that occurs for b ~ 1.2. On the other hand, model M2a appears to reproduce 
quite well the attractive short-distance part of PV^lb, b, b). The fact that model M2a gives a 
better estimate of PV^lb, b, b) than model Mia, while, at the same time, providing a slightly 
less accurate estimate of ^3 may seem contradictory at first sight. To understand it, let us 
note that (see Appendix for the derivation) 

xg{bi2)g{bi3)g{b23), (20) 

where the subscripts "pol" and " CG" refer to the polymer and the CG model, respectively, 
and g{b) is the zero-density center-of-mass pair distribution function, which is, by definition, 
identical in the polymer and in the CG model. Since g{b) is small for small values of b, 
keeping also into account that dhi2dbi^ gives a factor &i2^i3' small-distance behavior of 
V3(6i2, &13, ^23) is irrelevant for the computation of ^3. Hence it is much more interesting 
to compare 

F^{h) = 64(e-/3^3(W) _ 1)^3(6). (21) 

Such a quantity is reported in Fig. H] and shows that the relevant region corresponds to 
1 ~ ^ ^ 3. Moreover, while -^3(6) is mostly negative for polymers, we have -F3(&) = for 
model Mia, and, even worse, F^[b) > for model M2a. Hence, the compressible model gives 
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FIG. 4: Function -^3(6) as a function of 6, for polymers (FMa), for model M2a, and for the tetramer 
model of Ref. 22 (t). 

a correction to the results for model Mia which has the wrong sign: hence, the discrepancy 
with the polymer results increases. For comparison, we also include the estimate of Fj,{\)) 



for the tetramer model of Ref. |2^, which is quite close to the polymer result, confirming that 
the tetramer model is a very good CG model in the dilute regime. 

In Fig. [3] we also consider model M2b and the corresponding full monomer quantity 
(FMb). Potential FMb shows a clear logarithmic divergence as 6 — )■ 0, and indeed a general 
theoretical argument^>^ predicts V3(6, 6, 6) ~ —0.248 ln(l/6) for 6 — > 0. A fit of the small- 
distance data to V3 ~ —0.248 ln(a/6) gives a = 1.17. The corresponding expression gives 
a good fit of the full-monomer data for V3 up to 6 ^ 0.9. The three-body potential V3 
for model M2b is significantly different from the polymer one. For instance it is finite for 
b^O: 1^3(0, 0, 0) ^ -0.1 for model M2b. Clearly, the CG model M2b is unable to correctly 
reproduce the three-body interactions. 
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TABLE II: Compressibility factor Z(^) for the models introduced in Sec. IIIH for the tetramer 
model (t) of Ref. ,.22;. and for polymers (p) in the scaling limit. 

<!> p Mia Mlb M2a M2b M2c t 



0.135 1.187 1.18458(1) 1.17869(1) 1.18455(1) 1.18090(1) 1.18630(1 

0.27 1.393 1.38167(1) 1.36439(1) 1.38084(2) 1.36758(1) 1.38379(2 

0.54 1.854 1.80067(1) 1.74840(1) 1.79770(2) 1.75010(2) 1.80119(2 

0.81 2.371 2.23911(1) 2.14190(1) 2.23498(2) 2.13765(2) 2.23694(2 

1.09 2.959 2.70461(1) 2.55534(1) 2.70088(2) 2.54008(2) 2.69982(2 

2.18 5.634 4.55607(2) 4.18703(1) 4.56959(4) 4.08443(5) 4.55094(3 

4.36 12.23 8.29709(2) 7.47886(3) 8.36007(5) 6.9679(1) 8.31841(4 



1.18993(4) 
1.39852(6) 
1.8499(1) 



2.9090(1) 
5.2660(2) 
10.2056(1) 



B. The semidilute regime 

We now analyze the behavior of the models in the semidilute regime. For this purpose 
we have performed simulations of the CG models for several values of the polymer volume 
fraction 

<f = fR^c, (22) 

where c = L/V is the number of polymers per unit volume. We have determined the 
compressibility factor Z = f3P/c, the inverse compressibility K = d{cZ)/dc, the finite- 
density adimensional distribution P(cr, $) of Vg (a = Vg/Rg as before), and the ratio Sg{(^) = 
i?^($)/-R^. The compressibility factor was computed by using the molecular virial route^>^ 
and checked by comparing it with the (significantly less precise) result obtained by using 



the compressibility route (we compute K as described in Ref. |41| ). Both methods give the 
same results within errors, confirming our final estimates. For models Mia and Mlb we 
checked the Monte Carlo results using integral-equation methods — the hypernetted chain 
closure^ for model Mia and the Rogers- Young closure^ for model Mlb: we found very good 
agreement, indicating that these methods are very accurate for these soft-sphere models. 

In Table HTl and Fig.|5]we compare the estimates of Z for the different models. It is evident 
that considering radii dependent potentials is irrelevant for the thermodynamics, as already 
observed in the discussion of A^: models Mia, M2a, and M2c give completely equivalent 
estimates and so do models Mlb and M2b. Second, the CG model which uses the center of 
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FIG. 5: Compressibility factor Z as a function of ^. On the left we report results for models Mia, 
M2a, M2c, on the right we report the results for models Mlb, M2b. They are compared with the 
polymer prediction ZpM (full line, FM) (from Ref. 
lOOiZ/ZpM - 1). 
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I. In the insets we report the deviations 



mass as reference point appears to be more accurate than that using the central monomer. 
For $ = 1.09, which is the expected boundary of applicabihty of single-blob models, models 
Mia and M2a predict Z with an error of 9%, while model M2b underestimates Z by 14%. 

It is also interesting to compare the intermolecular distribution functions. In Fig. [6] 
we report the center-of-mass pair distribution function gciviib) and the polymer- midpoint 
pair distribution function gMp{b)- In the center-of-mass case (left panel) all results are in 
reasonable agreement for 6 > 1, while larger discrepancies are observed for 6 — > 0. For the 
midpoint distribution all models give similar curves, even for small values of b. This is due 
to the logarithmic divergence of the potentials, which enforces the condition gMp{b) — > for 
6 — )■ in all models. 

Finally, in Fig. [7| we report the distribution of the radius of gyration in the semidilute 
regime. For $ = 1.09 both types of coarse-graining reproduce correctly the polymer dis- 
tribution of Tg. For $ = 4.36 deviations are significantly larger. However, this should not 
be considered as a problem of the method, but rather as consequence of the single-blob 
model, which is only expected to work in the dilute regime. To perform a more quantitative 
comparison we consider the ratio 

S,{<^) = {rlh/Rl (23) 
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FIG. 6: Intermolecular distribution function for several models at ^> = 1.09 and 4.36 (the corre- 
sponding function is shifted upward for clarity). On the left we report results for models Mia, M2a, 
M2c, and the polymer center-of-mass distribution from full-monomer simulations (FMa); on the 
right we report the results for models Mlb, M2b, and the polymer distribution function associated 
with the polymer midpoint (FMb). 




FIG. 7: Distribution P{a,^) of cr = rg/Rg for $ = 1.09 and 4.36 for CG models M2a, M2b, and 
M2c and for polymers (FM). In the insets we report the deviations ASg = 100{Sg/ Sg^pM — 1), 
where Sg is the ratio (p3|) for the CG models and Sg^pM is the corresponding quantity for polymers. 

where the average is performed at polymer volume fraction $. For $ = 1.09 we obtain 
Sg{^) = 0.94161(5) (M2a), 0.92300(6) (M2b), 0.93630(6) (M2c), to be compared with the 
polymer result Sg{^) = 0.9238(2). In all cases differences are small, although model M2b ap- 
pears to reproduce better the polymer results. For larger values of $ discrepancies increase, 
see the insets in Fig. [71 but this is not surprising for a single-blob model. 
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IV. CONCLUSIONS 



In this paper we have performed a detailed study of single-blob models^i which are 
characterized by a fluctuating sphere radius and by density-independent potentials which 
are such to reproduce the radius-of-gyration distribution for an isolated polymer and the 
center-of-mass (or midpoint) polymer distribution function in the limit of zero polymer 
density. The results have been compared with those obtained in simpler single-blob models 
with fixed blob radius with the purpose of understanding if CG models based on compressible 
blobs provide a more accurate description of the thermodynamic behavior than fixed-radius 
CG models. Since we use single-blob models, this comparison can only be made in the dilute 
regime $ < 1, in which polymer overlaps are rare. 

As far as the thermodynamics and the intermolecular structure are concerned, we find 
that CG models that use the radii-dependent potentials behave no better than those in which 
polymers are represented by fixed-size spheres. However, the models with compressible soft 
spheres allow one to reproduce correctly (at least in the dilute regime in which single-blob 
models are expected to be predictive) the density dependence of the radius of gyration. For 
this quantity, the midpoint representation gives more accurate results. 

The fact that models with fluctuating sphere size and models with fixed sphere size 
have the same thermodynamic behavior is probably related to the fact that the radius of 
gyration is not the relevant length scale at finite density. If one wishes to improve the 
accuracy of the model, one should take into account overlaps which are characterized by a 
different length scale, the de Gennes-Pincus correlation length ^, which scales as ^ ~ ^g^~"' 
with 7 = z//(3z^ — 1) ^ 0.77 in the semidilute limit." To describe the semidilute regime, 
only multiblob approaches appear to be viable coarse-graining methods. In this respect, 
on the basis of the present results, we do not expect the approach of Vettorel et alM- 
to be more accurate for the thermodynamics than the more straightforward approach of 
Refs. 
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considering compressible soft blobs should not provide a model which gives a 
more accurate description of the polymer thermodynamics than those in which the blob size 
is fixed. Indeed, both representations equally fail to take into account blob overlaps. Of 
course, a compressible soft-blob model would reproduce better the density dependence of 
some structural properties, like the average radius of gyration and the form factor. 

One of the difficulties of multiblob approaches is the determination of the intra- and 
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inter-molecular potentials.— i^i^ For models with compressible blobs, a direct numerical 
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determination appears unfeasible, hence phenomenological approaches must be used. Ref. 
proposed a simple parametrization for the intermolecular blob potential, Eq. fllOp . For the 
single-blob case we find that this parametrization is quite accurate: model M2c is essentially 
equivalent to model M2a. 

Finally, we have compared the results for two different types of CG models. As for the 
thermodynamics, CG models for linear polymers that use the center of mass as reference 
point are significantly better than those that use the central monomer. Indeed, the estimates 
of the third virial coefficient and of the compressibility factor corresponding to models Mia 
and M2a are significantly closer to the full-monomer estimates than those obtained by using 
models Mlb and M2b. 
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Appendix A: Explicit expression for the third virial coefficient 

In this Appendix we wish to derive an expression for the third virial coefficient which 
explicitly depends on the three-body potential of mean force defined in Eq. f lTOj) . We consider 
a generic system of molecules with internal degrees of freedom and on each molecule we 
select a reference point S. Then, we indicate with (■)r the average over all internal degrees 
of freedom of the molecule such that the position of point 5* is r. The formalism can be 
applied both to polymers — in this case the average is over all polymer conformations and S 
can be taken as the polymer center of mass or the central monomer — and to the single-blob 
CG model — in this case S is the position of the sphere and the average is over all values of its 
radius. Analogously, we define (■)ri,r2 as the average over all additional degrees of freedom of 
two isolated molecules such that Si is in ri and 5*2 is in r2, and similar expressions involving 
three molecules. Finally, we define the zero-density S'-related pair distribution function 

gs{r) = (e-^^-)o,r, (Al) 
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where U12 is the intermolecular potential energy, and the corresponding correlation function 

hsir) = gsir) - 1 = (e"^^- - l)o,r = (/i2)o,r, (A2) 

where /12 is the usual Mayer function. 

The third virial coefficient for the model can be written as^ 

B, = ~h-Tu (A3) 

where 



',ri2,ri3; 



Tl = J dh^2d^rr3 ((/l2/l3)o,ri2,ri3 - (/l2)o,ri2(/l3)o,ri3) • (A4) 

Using definition ( !T9|) we can rewrite 

(/l2/l3/23)0,r,„r,3 = (g-'^^^^'^- '^"'^-^ - l) gs{ru) 9s{ri3) 9s{r23) 

+hs{ri2)hsirr3)hsir23) 

-((/12/13) o,ri2,ri3 — hs{r i2)hs{r 13) + 2 permutations) , (A5) 

where = ri3 — r23 and the two permutations correspond to replacing once 13 with 23 and 
the second time 12 with 23. Using this expression we end up with 

= -\j d'rud'r,, (e-/^^3(ri2,r,3,r23) _ 1) gsir,2)9sirn)9sir23) 

~ j dhi2dhi3hs{ri2)hs{ri3)hs{r23)- (A6) 

It is interesting to observe that this expression is identical to the one which applies to a 
system of monoatomic molecules (without additional degrees of freedom) interacting by 
means of a two-body and of a three-body potential. 

Using Eq. (1A6P we can compute the difference -B3,poi — B^ cg of the third virial coefficient 
for polymers and for the CG model. Since the S-related pair distribution function is, by 
definition, the same in the two models, the difference is given by Eq. f l20|) reported in the 
main text with g{r) = 9s{i^)- 
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